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PREFACE 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units 
verted  to  metric  (SI) 
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To  Obtain 
metres 
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CAVITY  DETECTION  AND  DELINEATION  RESEARCH 

ACOUSTIC  RESONANCE  AND  SELF-POTENTIAL  APPLICATIONS: 

MEDFORD  CAVE  AND  MANATEE  SPRINGS  SITES,  FLORIDA 

PART  I :  INTRODUCTION 

Background 

1.  One  of  the  more  challenging  problems  in  the  geophysical  field 
is  the  detection  and  delineation  of  subsurface  cavities.  The  presence 
of  either  naturally  developed  features  such  as  solution  cavities  in 
karstic  rock,  or  man-made  tunnels  in  rock  or  soil,  can  be  a  major  con¬ 
cern  in  a  variety  of  civil  and  military  applications.  Some  examples 
include  the  stability  of  critical  structures  such  as  dams  or  nuclear 
power  plants  founded  on  karst  materials,  the  erosion  of  foundation  mate¬ 
rials  by  piping,  and  clandestine  tunneling  beneath  or  near  military 
fortifications.  Most  major  construction  projects  include  a  program  of 
exploratory  borings,  some  of  which  may  detect  cavities  through  either 
accident  or  design.  The  high  cost  of  drilling  often  prohibits  the  use 
of  closely  spaced  borings  to  detect  and  delineate  subsurface  cavities. 
Ideally,  drilling  would  be  used  to  confirm  anomalies  detected  by  more 
economical  means,  possibly  geophysical  techniques. 

Purpose 

2.  The  purpose  of  this  study  was  to  assess  the  performance  of 
three  surface-applied  geophysical  methods  in  detecting  and  delineating 
known  subsurface  cavities.  Sites  with  known  dry  (air-filled)  or  wet 
(water-filled)  cavities  were  to  be  field-tested  to  determine  the  sensi¬ 
tivity  and  resolution  of  each  method  used. 
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Scope 

3.  This  report  documents  the  results  of  acoustic,  sonar,  and 
self-potential  (SP)  tests  carried  out  at  two  karst  sites  in  Florida. 
Medford  Cave,  a  dry  site,  was  tested  using  the  acoustic  method.  Manatee 
Springs,  a  wet  site,  was  tested  using  sonar  and  SP  methods.  Locations 
of  these  sites  are  shown  in  Figure  1.  Field  test  procedures,  the  data 
obtained,  and  an  analysis  of  results  for  both  sites  are  presented  herein. 


5 


PART  II :  METHODOLOGY 


Acoustic  Resonance 


4.  Acoustic  resonance  techniques  for  soil  studies  were  first 
employed  by  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES) 

in  the  period  1975-76  (Ballard,  1977).  This  work  was  directed  to  identi¬ 
fying  shallow  man-made  tunnels  in  soil,  for  military  purposes.  It  was 
believed  that  the  course  of  a  shallow  tunnel  in  soil  could  be  traced  by 
measuring  the  ground  surface  wave  amplitudes  produced  by  a  loudspeaker 
placed  in  the  tunnel  mouth.  In  this  procedure,  the  speaker  was  first 
tuned  to  the  resonant  frequency  of  the  cavity,  and  wave  amplitudes  at 
preselected  locations  on  the  ground  surface  were  measured  with  a  very 
sensitive  geophone  probe.  A  schematic  of  the  typical  equipment  configu¬ 
ration  used  for  an  acoustic  survey  is  shown  in  Figure  2. 

5.  The  acoustic  resonance  technique  met  with  success  in  the 
initial  testing  over  lateral  distances  of  a  few  hundred  feet  in  soil, 
and  was  later  used  to  investigate  shallow-buried  natural  solution  chan¬ 
nels  in  rock  (Cooper  and  Bieganousky,  1978).  For  natural  or  man-made 
features  at  shallow  depths,  the  survey  is  usually  carried  out  by  measur¬ 
ing  induced  wave  amplitudes  at  the  ground  surface  at  points  in  a  pre¬ 
determined  grid  pattern.  Grid  increments  equal  to  about  one-half  the 
estimated  diameter  of  the  features  of  interest  are  typically  employed. 

The  general  rationale  used  for  interpretation  of  results  is  that  the 
maximum  amplitude  of  surface  waves  emanating  from  shallow  cavities  should 
occur  immediately  above  the  cavity  roof,  and  that  the  absolute  magnitude 
of  the  induced  signal  will  diminish  (through  attenuation)  with  increas¬ 
ing  distance  from  the  loudspeaker  location. 

6.  To  date,  this  technique  has  been  applied  only  at  locations 
where  the  subsurface  cavity  was  accessible  through  an  opening  at  the 
ground  surface.  However,  it  would  also  be  possible  to  use  a  high- 
intensity  borehole  acoustic  source  to  excite  cavities  inaccessible  by 
other  means.  The  most  recent  application  of  acoustic  resonance  methods 
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was  at  the  Medford  Cave  test  site,  and  results  of  this  study  are  pre¬ 
sented  later  in  this  report. 


Sonar 


7.  The  acoustic  resonance  technique  previously  described  may 
only  be  used  in  air-filled  cavities.  For  water-filled  cavities,  the 
WES  desired  to  use  a  similar  approach  by  substituting  a  sonar  source  for 
the  acoustic  loudspeaker.  High- intensity  sonar  sources  suitable  for  the 
purpose  have  been  developed  by  the  U.  S.  Naval  Research  Laboratory  (NRL) 
at  Orlando,  Fla.  NRL  sonar  sources  have  been  extensively  tested  by  the 
U.  S.  Navy,  other  Government  agencies  and  private  firms,  and  are  docu¬ 
mented  in  the  literature  (Groves,  1974).  These  devices  are  packaged  in  a 
variety  of  sizes  to  suit  specific  applications;  however,  one  of  the 
limiting  factors  is  depth  of  submergence.  This  limitation  did  prove  to 
be  a  problem,  as  will  be  discussed. 

8  It  was  desired  to  test  the  feasibility  of  using  sonar  sources 
to  detect  water-filled  cavities  located  at  depths  up  to  100  ft*  below 
the  ground  surface.  Surface  delineation  of  the  cavity  feature(s)  would 
be  accomplished  using  the  same  techniques  employed  in  the  acoustic 
resonance  study  of  air-filled  cavities.  One  of  the  more  important  ques¬ 
tions  to  be  answered  was  whether  the  sonar  source  could  deliver  suffi¬ 
cient  energy  to  propagate  measurable  signals  for  distances  of  at  least 
several  hundred  feet  in  soil  and  rock.  A  feasibility  study  using  a 
sonar  source  has  been  conducted  at  the  Manatee  Springs  test  site,  and 
results  are  presented  later  in  this  report. 

Self-Potential 


9.  The  movement  of  pore  water  in  subsurface  strata  and  spatial 
variations  in  the  chemical  composition  of  the  strata  and  pore  fluid  can 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  3. 
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result  in  the  production  of  relatively  low-level  direct  current  (DC) 
earth  currents,  known  as  self-potentials  (SP) .  The  natural  steady-state 
electrical  field  observed  at  the  ground  surface  may  range  in  magnitude 
from  a  few  millivolts  to  several  hundred  millivolts  or  more,  depending 
on  local  conditions  and  the  method  of  measurement.  For  more  than  50 
years,  the  surface  distribution  of  the  SP  field  has  been  used  to  locate 
ore  bodies,  and  more  recent  applications  include  identification  of 
seepage  areas  in  dams  and  reservoirs,  and  geothermal  studies  (Corwin 
and  Hoover,  1979,  and  Ogilvy,  Ayed ,  and  Bogoslovsky,  1969). 

10.  For  the  case  of  nongeothermal  water  flowing  through  under¬ 
ground  channels,  the  flowing  fluid  forms  a  double-boundary  layer  at  the 
channel  wall  which  results  in  a  net  negative  electrical  charge  in  the 
surrounding  mineral  crystals  as  shown  in  Figure  3  (Corwin  and  Hoover, 
1979) .  If  the  flowing  water  is  fresh  (has  relatively  low  salinity) , 
the  negative  potential  may  be  of  sufficient  magnitude  to  produce  a 
significant  variation  in  the  surface  SP  field. 

11.  The  usual  procedure  for  determining  SP  distribution  at  the 
surface  is  to  first  place  a  reference  electrode  at  a  relatively  large 
distance  from  the  survey  area  (usually  upstream,  if  the  direction  of 
fluid  flow  is  known  with  reasonable  confidence).  A  second  (survey) 
electrode  is  then  traversed  along  a  predetermined  grid  pattern,  making 
repeated  DC  potential  measurements  between  the  reference  and  survey 
electrodes.  Measurements  can  be  made  using  a  voltmeter  with  high 
internal  resistance  (say,  50  megohms  or  more)  so  as  to  minimize  current 
flow  in  the  circuit.  This  is  necessary  to  prevent  undesirable  polari¬ 
zation  of  the  measurement  electrodes.  Sufficient  data  are  obtained  to 
adequately  define  the  surface  SP  field,  the  distribution  of  which  may 
indicate  subsurface  cavities.  An  investigation  of  this  type,  conducted 
at  the  Manatee  Springs  test  site,  is  described  herein. 
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PART  III:  DATA  ACQUISITION 


Medford  Cave  Test  Site  (Dry  Site) 

12.  The  Medford  Cave  test  site  is  located  in  gently  rolling, 
partly  wooded  terrain,  which  has  local  relief  in  the  site  vicinity  of 
about  100  ft.  The  site  area  is  characterized  by  a  thin  (typically  1  to 
5  ft  thick)  layer  of  sand  overlying  limestone. 

13.  Solution  activity  in  the  limestones  of  the  area  has  produced 
karst  topography  typified  by  sinkholes,  subsurface  cavities,  and  little 
or  no  surface  drainage.  The  top  of  rock  is  typically  pinnacled,  with 
filled  clay  pockets  between  pinnacles.  Cap  rock  at  the  site  is  a  3-ft- 
thick  layer  of  hard  molluscan  limestone  of  the  Hawthorne  formation. 

This  cap  is  unconformably  underlain  by  limestone  of  the  Ocala  formation, 
which  is  soft  to  very  soft,  friable,  and  solutioned  to  form  the  known 
portions  of  the  Medford  Cave  system.  The  water  table  at  the  site  is  at 
a  depth  of  about  115  ft  below  the  ground  surface. 

14.  A  plan  view  of  the  site,  showing  the  general  outline  of  the 
complex  cavity  system,  profile  lines,  and  borings  is  presented  in 
Figure  4.  Typical  cross  sections  through  the  main  cavity  system  are 
shown  in  Figure  5.  A  subsurface  profile  along  the  80  north  grid  line  is 
shown  in  Figure  6.  The  cavity  outline  shown  in  Figure  4  is  only  approxi¬ 
mate  since  the  complexity  of  the  features  precludes  complete  mapping; 
outlines  of  uncertain  accuracy  are  indicated  by  dashed  rather  than  con¬ 
tinuous  lines.  The  original  cavity  survey  was  performed  by  Southwest 
Research  Institute  (SwRI)  in  conjunction  with  the  Florida  Department  of 
Transportation.  Minor  corrections  to  the  plotted  outlines  were  later 
made  by  WES. 

15.  It  can  be  seen  from  Figures  4  and  5  that  the  cavity  system 
is  typically  encountered  at  depths  ranging  from  10  to  50  ft  below  the 
ground  surface,  and  that  the  cavities  are  quite  irregular  in  shape, 
cross  sectional  area,  and  principal  direction.  This  degree  of  complexity, 
although  probably  characteristic  of  most  solution-developed  features,  is 
clearly  a  severe  test  of  resolution  for  any  detection  and  delineation 


method  employed.  Even  so,  sufficient  control  at  the  site  exists  to  make 
a  meaningful  determination  of  resolution  for  a  range  of  feature  sizes. 

16.  Considering  that  features  from  a  few  feet  to  30  ft  or  more 
in  width  were  to  be  detected,  a  surface  measurement  increment  of 

2.5  ft  was  adopted  as  a  reasonable  choice  for  the  acoustic  survey. 

Three  areas  selected  for  survey  were  known  to  overlie  small,  interme¬ 
diate,  and  large  cavity  features  and  are  identified  as  areas  1,  2,  and 
3,  respectively,  in  Figure  7.  Areas  2  and  3  overlap,  and  were  surveyed 
using  different  locations  of  the  loudspeaker  within  the  cavity  system. 
Access  to  the  speaker  locations  was  through  the  main  cavity  entrance,  in 
the  case  of  areas  2  and  3,  and  a  secondary  cavity  opening  in  area  1. 

The  loudspeaker  used  was  a  commercially  available  unit  with  a  12-in. - 
diam  bass  speaker  capable  of  handling  60  watts  root-mean-square  (RMS) 
continuous  sine  wave  electrical  input  signals.  The  speaker  was  excited 
with  steady-state  sine  waves  using  a  signal  generator  and  a  30-watt  RMS 
rated  high  fidelity  amplifier.  After  the  speaker  was  placed  at  the 
desired  location,  a  frequency  sweep  was  made  to  determine  the  acoustic 
frequency  which  produced  maximum  amplitude  ground  motion  (apparent 
resonance).  In  the  case  of  survey  areas  1,  2,  and  3,  the  apparent 
resonant  frequencies  were  116,  108,  and  102  Hz,  respectively.  During 
each  survey  the  initial  frequency  and  power  settings  were  held  constant. 

17.  As  shown  in  Figure  4,  a  primary  spacing  of  20  ft  was  used  in 
laying  out  the  survey  grid.  The  bulk  of  the  acoustic  data  was  taken  at 
2.5-ft  intervals  along  the  primary  grid  lines,  but  some  diagonal  line 
and  feature  tracing  measurements  were  also  made  where  site  conditions 
permitted  (parts  of  the  northern  half  of  the  survey  area  are  rather 
heavily  wooded,  which  complicated  the  measurement  process) .  Acoustic 
measurements  were  made  with  a  portable  geophone  system  consisting  of 
three  series  wired  vertical  coil  elements  rigidly  housed  in  a  hand- 
carried  spiked  probe.  A  portable  battery-powered  amplifier  and  oscillo¬ 
scope  were  used  to  measure  the  peak-to-peak  voltage  of  the  amplified 
signal  from  the  probe.  When  surveying,  the  probe  was  pushed  into  the 
soil  until  the  spike  was  firmly  embedded  and  then  the  signal  magnitude 
was  read  from  the  oscilloscope  with  the  operator  standing  still  and  not 


touching  the  probe.  Amplifier  gain  settings  were  held  constant  through¬ 
out  the  surveys  (nominally  at  1.5  times  the  geophone  output).  Approxi¬ 
mately  20  sec  was  required  to  make  a  measurement  at  each  station  location. 
A  contour  plot  of  acoustic  signal  magnitude  for  test  area  1,  the  loud¬ 
speaker  location,  and  the  estimated  (SwRI)  and  verified  (WES)  cavity 
location  in  this  locale  are  shown  in  Figure  8.  A  tabulation  of  the 
magnitude  of  the  acoustic  signal  detected  at  each  measurement  location  is 
presented  in  Table  1.  The  measured  signal  magnitudes  are  expressed  in 
peak-to-peak  volts  at  the  resonant  frequency  of  excitation.  Since  the 
intent  of  this  survey  was  to  measure  relative  rather  than  absolute  signal 
magnitudes,  no  rigorous  calibration  of  the  surface  probe  geophones  was 
carried  out.  However,  a  sensitivity  value  of  10.0  volts/in. /sec  is 
approximately  (+30  percent)  correct  for  the  gain  setting  used. 

Figures  9  and  10  show  contour  plots  of  signal  magnitude  and  loudspeaker 
locations  for  survey  areas  2  and  3,  respectively.  The  acoustic  data  for 
area  2  is  tabulated  in  Table  2  and  for  area  3  in  Table  3.  All  of  the 
acoustic  data  were  obtained  by  a  two-man  crew  in  the  period  18-20  August 
1979. 


Manatee  Springs  Test  Site  (Wet  Site) 

18.  The  Manatee  Springs  test  site  is  located  approximately  6 
miles  west  of  Chiefland,  Fla.,  on  the  east  bank  of  the  Suwanee  River, 
as  shown  in  Figure  1.  The  site  is  within  the  Pamlico  marine  terrace 
which,  in  the  site  area,  consists  of  a  broad  shelf  of  Tertiarv  limestones 
covered  with  a  thin  veneer  of  Pleistocene  sand.  Elevation  in  the 
vicinity  of  the  site  is  10  to  15  ft,  and  relief  on  site  is  on  the  order 
of  3  ft.  Erosion  of  the  limestone  prior  to  deposition  of  the  sand  had 
produced  an  irregular  limestone  surface,  as  evidenced  by  a  13-ft 
variation  in  elevation  of  the  top  of  rock  (from  borings  on  site).  The 
coastal  terrace  is  characterized  by  karst  topography  which  includes 
sinkholes,  vertical  piping,  poor  surface  drainage,  and  solution-developed 
cavity  systems.  Large  solution  features  on  site  are  encountered  in  the 
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Ocala  and  Wllliston  limestone  units  of  Eocene  age,  at  depths  ranging 
from  45  ft  (spring  boil)  to  90  ft  (main  cavity  system). 

19.  A  plan  view  of  the  Manatee  Springs  site  is  shown  in 

Figure  11.  In  Figure  12,  the  general  outline  of  the  WES-mappped  portion 
of  the  cavity  is  shown;  dashed  outlines  indicate  areas  not  mapped  for 
this  study.  Also  shown  in  Figure  12  is  the  survey  grid  used  to  make 
surface  measurements  and  the  locations  of  exploratory  borings.  In  the 
initial  phase  of  developing  the  test  site,  it  was  necessary  to  locate 
the  subsurface  features  with  greater  precision  than  that  afforded  by 
existing  maps  (Figure  11).  Borings  S-l  and  C-l  were  drilled  to  provide 
reference  points  from  which  more  precise  measurements  of  the  cavity 
features  could  be  made.  Members  of  the  Cave  Divers  Section,  National 
Speleogical  Society  then  resurveyed  the  cavity  features  of  interest  so 
that  other  boreholes  (C-2  through  C-5)  could  be  situated  to  bracket  a 
side  tunnel  feature  for  purposes  of  a  different  study  not  reported 
herein.  Boreholes  E-l  through  E-6  and  UL-2  were  drilled  later  to  better 
define  areas  of  interest  at  the  site.  From  these  data  and  divers' 
observations,  the  cross  sections  shown  in  Figures  13  and  14  were 
developed.  It  should  be  stressed  that  the  main  cavity  outlines  shown  in 
Figure  12  are  correct  insofar  as  major  features  are  concerned,  but  that 
the  divers  reported  areas  where  secondary  cavity  features  are  so  numer¬ 
ous  that  the  limestone  resembles  Swiss  cheese.  However,  as  In  the  case 
of  Medford  Cave,  major  accessible  features  have  been  mapped  in  sufficient 
detail  to  premit  a  meaningful  test  of  resolution  for  the  survey  methods 
employed. 

20.  Sonar  and  SP  investigations  at  the  site  were  completed  in 
the  period  9-15  October  1980  by  a  three-man  WES  field  crew.  The  SP 
survey  was  conducted  by  taking  measurements  with  the  survey  electrode 
at  10-ft  intervals  along  the  E-W  grid  lines  already  established.  The 
reference  electrode  was  located  in  Sue's  Spring,  approximately  750  ft 
upstream,  as  shown  in  Figure  11.  The  potential  difference  between 
reference  and  survey  electrodes  was  measured  with  a  commercially  avail¬ 
able  digital  volt-ohmmeter  having  an  Internal  impedance  of  1000  megohms. 
The  electrodes  used  were  porous  pots  of  the  copper-copper  sulfate  type, 
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which  were  connected  to  the  volt-ohmraeter  with  multistrand  16  AWG  insu¬ 
lated  wire.  All  of  the  measured  SP  data  are  tabulated  in  Table  4. 

Figure  15  is  a  plan  view  of  the  SP  survey  area  which  shows  the  survey 
grid,  subsurface  cavity  outline,  and  equlpotential  contour  lines 
developed  on  20-mv  intervals. 

21.  In  the  sonar  feasibility  study,  the  sonar  source  was  posi¬ 
tioned  as  shown  in  Figure  16  using  a  small  boat  to  suspend  the  source  in 
the  sinkhole  identified  in  Figure  11  as  the  "Catfish  Hotel."  This  sink¬ 
hole  communicates  with  the  main  cavity  system  at  a  depth  of  about  90  ft 
below  the  ground  surface.  Because  of  previously  mentioned  limitations 
on  its  depth  of  submergence  (maximum  depth  =  75  ft),  the  sonar  source 
could  not  be  placed  directly  in  the  main  cavity.  At  depths  greater  than 
75  ft  the  internal  pressure-equalizing  system  ir.  the  device  cannot 
compensate  for  hydrostatic  pressure  on  the  moving  diaphragm,  causing  a 
drastic  decrease  in  output. 

22.  A  frequency  sweep  was  first  performed  to  determine  the  dis¬ 
crete  frequency  at  which  surface-received  signals  had  maximum  amplitude. 
The  sonar  source  was  driven  using  a  tone  burst  signal  generator  and  an 
amplifier  capable  of  delivering  250  watts  of  continuous  power.  However, 
the  sonar  source  could  accept  only  50  watts  power  input  so  this  excita¬ 
tion  level  was  maintained  throughout  testing.  After  determining  that 
123  Hz  was  the  best  signal  propagation  frequency  for  this  site,  the 
pulse  generator  was  used  to  automatically  pulse  the  sonar  source  (1  sec 
on,  1  sec  off)  during  data  acquisition.  Since  there  is  always  some 
background  seismic  noise  at  any  test  site,  pulsing  the  sonar  signal  was 
a  convenient  way  to  differentiate  between  noise  and  the  signal  of 
interest  when  weak  signals  were  being  received.  The  peak-to-peak 
amplitude  of  surface-arriving  sonar  pulse  was  measured  using  the  portable 
geophone  system  described  in  paragraph  17;  however,  the  amplifier  gain 
was  increased  to  yield  about  400  volts/in. /sec  total  system  output 
sensitivity  in  an  effort  to  detect  very  feeble  signals  at  some  locations 
on  site.  It  had  been  hoped  that  naturally  occurring  water  turbelence 

in  the  main  flow  channel  might  produce  measurable  signals  at  the  surface; 
however,  background  noise  surveys  indicated  that  these  signals,  if 
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present,  were  too  feeble  to  detect  even  with  the  increased  gain  settings 
used.  Hence,  only  sonar  pulse  data  are  presented  herein.  The  sonar 
data  are  tabulated  in  Table  5.  Results  of  the  sonar  feasibility  survey 
are  also  shown  in  Figure  17,  which  is  a  contour  plot  of  peak-to-peak 
voltage  amplitude  of  received  signals  using  1-volt  contour  intervals. 
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PART  IV:  ANALYSIS  OF  RESULTS 


23.  Both  the  acoustic  resonance  and  sonar  methods  have  to  date 
been  employed  as  a  means  to  delineate  cavities  whose  presence  is  known. 
Consequently,  an  evaluation  of  these  methods  is  properly  directed  to 
their  ability  to  delineate  cavity  features.  On  the  other  hand,  the 
self-potential  method  may  be  employed  to  both  detect  and  delineate 
subsurface  cavities,  provided  they  are  water-filled,  some  flow  is  tak¬ 
ing  place,  and  other  necessary  conditions  (such  as  low  salinity  pore 
fluid)  obtained.  SP  results  in  the  current  study  are  evaluated 
accordingly. 


Acoustic  Resonance  Survey  (Medford  Cave) 

24.  The  rationale  for  interpreting  acoustic  resonance  data 
from  simple  features,  as  for  example  a  shallow  tunnel  in  soil,  is  that 
the  maximum  amplitude  of  the  speaker- induced  wave  occurs  immediately 
above  the  cavity  feature.  This  has  proven  to  be  the  case  in  previous 
investigations,  but  the  features  studied  had  usually  been  situated  a 
few  feet  below  the  surface  and  the  maximum  horizontal  feature  dimen¬ 
sions  had  been  on  the  order  of  1  to  3  ft,  typically.  Medford  Cave 
generally  represented  a  much  more  complex  condition,  both  in  terms  of 
depth  of  burial  of  the  features  and  in  terms  of  their  size.  However,  in 
acoustic  survey  area  1,  the  cavity  environment  most  closely  approximated 
conditions  previously  studied  and  the  acoustic  amplitude  contour  plot 
shown  in  Figure  8  is  centered  on  the  known  feature  and  does  indicate  its 
general  extent.  It  had  been  recognized  that  larger  and  more  complex 
features  might  not  prove  to  be  susceptible  to  such  a  simplistic  approach, 
because  the  drastic  changes  in  geometry  could  easily  produce  wave  inter¬ 
ference  patterns,  signal  cancellation,  beat  frequencies,  etc.  It  was 
argued  that  the  signal  received  at  the  surface  should  attenuate  as  a 
function  of  range  from  the  source,  in  a  homogeneous  medium,  so  that 
departures  from  this  classic  spatial  pattern  might  be  used  as  a  means  of 
detecting  anomalous  features  (cavities). 
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25.  An  approach  of  this  type  was  tried  using  the  data  obtained 

in  acoustic  survey  area  1.  The  data  as  recorded  at  the  site  were  of  the 
form  (x,y,a)  where  x  and  y  represent  the  rectangular  grid  location 
of  the  measuring  instrument  and  a  the  amplitude  of  the  signal  recorded 
at  the  location.  The  data  were  transformed  to  the  form  (d,a)  where  d 
represents  the  surface  distance  from  the  point  directly  above  the  signal 
source  (106,260)  to  the  recording  instrument  location.  For  this  data 
set,  the  distances  ranged  from  1  to  61.61  ft,  and  the  signal  amplitude 
measured  peak-to-peak  ranged  from  0  to  132  mv. 

26.  A  least-squares  fit  was  made  to  the  attenuation  data  using 

various  functional  relations.  The  relation  a  =  A  +  B/d  was  selected 
as  the  best  on  the  basis  of  an  R-squared  (square  of  the  coefficient  of 
correlation  R)  value  of  0.41  and  the  fact  that  it  gave  reasonable 
predicted  values  over  the  range  of  data.  The  data  are  plotted  with 

the  fitted  curve  shown  in  Figure  18. 

27.  The  residuals  from  the  relation  a  =  A  +  B/d  are  not 
normally  distributed;  therefore,  Chebyshev's  inequality  (Hogg  and  Craig, 
1970),  Pr(|a  -  a|  _>  KS)  <  1/K  is  used  to  calculate  the  93-percent 
(corresponding  to  K  =  3  )  confidence  interval.  Chebyshev's  inequality 
states  that  the  random  variable  a  is  distributed  about  its  mean  value 
a-  according  to  the  above  relation  for  any  positive  number  K  ,  where  S 
represents  the  standard  deviation.  The  upper  and  lower  bounds  of  the 
confidence  interval  in  Figure  18  are  curves  fitted  to  the  points  calcu¬ 
lated  using  the  formula  from  Draper  and  Smith  (1966): 

d.  +  K  •  S  • 

l  — 

where 

dt  =  the  value  of  a  in  the  relation  a  =  A  +  B/d  corres¬ 
ponding  to  d^ 

K  =  3 

S  =  the  standard  deviation  of  sample  data 

T 

=  the  vector  (1,  d^) 
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X^  =  X^  transpose 


The  matrix  C  =  (X^X)-1 
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and  n  equals  the  number  of  data  points.  Upper-  and  lower-bound  curves 
for  93-percent  confidence  interval  are  also  shown  in  Figure  18.  A  total 
of  40  data  points  fall  above  the  upper-bound  curve,  and  these  data  were 
normalized  by  dividing  by  the  upper-bound  curve  value  at  the  same  range. 
The  resulting  set  of  40  data  points,  ranging  in  amplitude  ratio  from 
1.25  to  3.74,  identify  higher  signal  amplitudes  than  might  be  predicted 
from  the  upper-bound  confidence  curve.  Since  the  rationale  was  that 
anomalously  high  amplitude  signals  should  occur  over  a  cavity  feature(s), 
these  data  were  used  to  develop  the  contour  plot  shown  in  Figure  19. 
Comparing  Figures  8  and  19,  it  can  be  seen  that  the  plot  in  Figure  19 
appears  to  conform  more  closely  to  the  WES-verif ied  cavity  location  in 
acoustic  survey  area  1.  Figures  20  and  21  present  similarly  derived 
plots  of  acoustic  signal  attenuation  versus  range  from  the  source  for 
acoustic  survey  areas  2  and  3,  respectively.  Comparing  Figures  18,  20, 
and  21,  it  can  be  seen  that  considerably  more  data  scatter  is  evident  in 
Figures  20  and  21  than  in  Figure  18.  Hence,  it  might  be  expected  that 
the  more  complex  features  investigated  in  survey  areas  2  and  3  might  not 
be  imaged  so  well  by  this  method.  Figures  22  and  23  present  normal¬ 
ized  signal  attenuation  contour  plots  for  survey  areas  2  and  3,  respec¬ 
tively,  and  little  improvement  can  be  seen  when  comparing  these  with  the 
original  raw  data  contour  plots  (Figure  9  versus  22,  Figure  10  versus 
23).  It  appears  that  the  signal  attenuation  versus  range  approach  does 
not  offer  significant  advantages  when  cavity  features  are  both  large  and 
complex  and  are  situated  a  few  tens  of  feet  below  the  ground  surface. 

In  an  earlier  investigation  (Cooper  and  Bieganousky,  1978),  the  acoustic 
resonance  technique  had  been  successfully  used  to  trace  shallow  solution 
cavities  in  limestone.  Some  of  these  cavity  features  were  situated  at 
depths  comparable  to  the  cavities  at  Medford  Cave;  however,  the  Medford 
Cave  main  cavity  features  are  at  least  an  order  of  magnitude  larger  in 
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lateral  dimensions  than  the  known  cavity  features  previously  investigated. 
So,  it  does  appear  that  this  acoustic  resonance  technique  suffers  a 
marked  decrease  in  resolution  when  applied  to  large  and  complex  cavern¬ 
ous  features. 


Sonar  Study  (Manatee  Springs) 

28.  By  inspection  of  Figure  17,  it  can  be  seen  that  easily 
measurable  sonar  signals  could  be  propagated  to  distances  of  200  ft  or 
more  from  the  source.  This  satisfied  the  primary  question  as  to  the 
feasibility  of  the  approach.  However,  it  can  also  be  seen  from 
Figure  17  that  numerous  areas  of  high-amplitude  sonar  signals  were 
detected  at  points  remote  from  the  known  main  cavity  channel.  This 
result  was  unexpected,  and  unfortunately  time  and  cost  limitations  pre¬ 
cluded  drilling  these  locations  in  an  effort  to  develop  a  rationale  for 
cause  and  effect.  No  signal  attenuation  methodology  was  applied  since 
there  appeared  to  be  no  consistent  relationship  between  signal  amplitude 
and  cavity  location. 

29.  It  can  be  seen  in  Figure  17  that  the  right  half  of  the 
sonar  amplitude  contour  plot  is  comprised  of  very  low  level  values.  It 
is  believed  that  the  sonar  signals  did  ^ot  effectively  excite  the  main 
cavity  system  since  the  sonar  source  had  to  be  located  in  the  adjacent 
sinkhole  because  of  the  previously  discussed  limitations  on  submergence. 

Self-Potential  Survey  (Manatee  Springs) 

30.  In  a  situation  where  the  direction  of  water  flow  is  known, 
as  at  Manatee  Springs,  the  reference  electrode  may  be  positioned  up¬ 
stream  in  which  case  the  presence  of  a  flow  channel  should  certainly  be 
indicated  as  a  negative  anomaly  in  the  surface  SP  field.  The  geology  of 
the  Manatee  Springs  site  suggests  that  the  limestone  is  quite  pervious, 
and  in  fact  this  has  been  verified  by  visual  observations  by  cave  divers, 
so  that  a  negative  bias  in  the  SP  field  might  be  expected  from  water 
flowing  through  numerous  small  interconnected  passages  as  well  as  in 
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major  solution  features  within  the  survey  zone.  In  fact,  the  raw  SP 
contour  plot  presented  in  Figure  15  shows  only  one  localized  area  where 
pronounced  positive  SP  values  were  recorded. 

31.  In  an  attempt  to  gain  better  resolution  in  the  SP  survey 
plot,  it  was  decided  that  the  mean  regional  SP  might  be  removed  by  sub¬ 
tracting  the  average  SP  value  from  each  SP  measurement.  The  average  SP 
value  was  determined  to  be  -25  mv,  and  using  this  value,  a  revised  SP 
contour  plot  was  developed  as  shown  in  Figure  24.  Comparing  Figure  15 
(raw  SP)  with  Figure  24  (residual  SP) ,  it  can  be  seen  that  Figure  24 
provides  a  spatial  representation  which  is  simpler  and  conforms  better  to 
the  known  location  of  cavity  features.  The  zero  SP  residual  contour 
area  appears  to  crudely  delineate  the  major  zone  of  flow,  although  an 
exact  delineation  is  probably  not  to  be  expected  when  the  target  system 
is  located  90  ft  below  the  ground  surface.  Only  the  areal  extent  of  the 
system  can  be  predicted  using  SP  techniques  developed  to  date;  the  pre¬ 
diction  of  depth  is  not  yet  within  the  state  of  the  art. 

Summary  and  Conclusions 

32.  It  is  concluded  that  the  acoustic  resonance  technique  is 
useful  in  delineating  shallow  (depth  less  than  20  ft)  subsurface  cavity 
features  whose  maximum  horizontal  dimensions  are  less  than  10  feet,  and 
that  the  interpretation  in  such  cases  is  reasonably  straightforward. 
However,  as  it  usually  the  case  with  most  geophysical  techniques,  in¬ 
creasing  geologic  complexity  and  depth  of  burial  have  drastic  adverse 
effects  on  resolution.  A  sonar  feasibility  study  was  also  conducted  and 
it  was  demonstrated  tiiat  a  sonar  source  could  deliver  adequate  energy 

to  excite  a  deep  water-filled  cavity  system;  however,  the  known  cavity 
features  at  Manatee  Springs  could  not  be  delineated  using  the  limited 
data  obtained.  The  only  NRL  sonar  source  available  at  the  time  of  this 
study  was  not  well  suited  to  site  conditions,  so  this  method  may  deserve 
further  attention. 


33.  In  this  study  the  SP  method  was  applied  to  the  problem  of 
detection  and  delineation  of  water-filled  solution  cavities  in  which 


flow  was  occurring.  A  known  large  solution  channel  located  at  a  depth 
of  90  ft  was  detected  by  surface  methods,  and  a  crude  delineation  of  the 
feature  was  also  achieved.  There  is  as  yet  no  technique  for  determining 
feature  depth  using  surface  SP  methods.  Improvement  in  cavity  feature 
resolution  was  obtained  in  one  phase  of  the  Medford  Cave  acoustical 
survey  and  in  the  Manatee  Springs  SP  survey  by  detrending  the  raw  data 
based  on  (a)  hypotheses  as  to  the  physical  processes  involved,  and 
(b)  careful  observation  of  the  apparent  trends  in  the  raw  data.  However, 
these  techniques  may  not  be  applicable  under  all  conditions. 
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NOTE:  Measured  signal  amplitudes  are  expressed  in  peak-to-peak  volts  read  from  an  oscillo¬ 
scope.  To  convert  these  readings  to  approximate  peak-to-peak  particle  velocity  in 
inches  per  second,  divide  the  voltage  value  by  10  (see  text,  page  11). 
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Se If -Pot  ential  Survey  Data.  Manatee  Springs 


NOTE:  Readings  are  expressed  in  millivolts  DC  read  between  the  survey  and  reference  electrodes. 


Sonar  Data.  Manatee  Springs  Site 


NOTE:  Measured  signal  amplitudes  are  expressed  in  peak-to-peak  volts  read  from  an  oscilloscope.  To  convert 
these  readings  to  approximate  peak-to-peak  particle  velocity,  divide  the  voltage  value  by  400  (see 
text,  page  13). 
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Figure  4.  Plan  view  of  the  Medford  Cave  test  site  showing  cavity 
system,  boring,  and  profile  locations 
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Figure  5.  Typical  cross  sections  through  the  Medford  Cave 

main  cavity  system 
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Figure  6.  Subsurface  profile  along  the  80  north  grid  line 


300,280 


300,0 


Figure  7.  Location  of  acoustic  resonance  survey  areas  with  respect  to 
the  main  cavity  system,  Medford  Cave 


Figure  9.  Contour  plot  of  acoustic  signal  magnitude  showing  location 
of  acoustic  source,  test  area  2,  Medford  Cave 
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Figure  10.  Contour  plot  of  acoustic  signal  magnitude  showing 
location  of  acoustic  source,  test  area  3,  Medford  Cave 
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Figure  12.  Location  of  WES  survey  area  at  Manatee  Springs  test  site  showing 
boring  locations  and  general  outline  of  main  cavity 
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Figure  13.  Cross  section  between  borings  C-l  and  S-l 
showing  main  cavity  system 
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Figure  14.  Cross  section  through  borings  C-2  to  C-5 
showing  secondary  cavity 
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Figure  16.  Cross  section  showing  sonar  source  location  in  sinkhole 
conununicat ing  with  main  cavity  system 
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Acoustic  signal  attenuation  versus  range  from  the  source 
survey  area  1,  Medford  Cave 
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Figure  21.  Acoustic  signal  attenuation  versus  range  plot  for  test 
area  3,  Medford  Cave  test  site 


Figure  22.  At tenuat ion-no rma ..ized  acoustic  contour  plot  for  test 
area  2,  Medford  Cave  test  site 
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Figure  23.  Attenuation  normalized  acoustic  contour  plot  for 
test  area  3,  Medford  Cave  test  site 
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